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Abstract 

> 

Three variational approaches, the hyperspherical-harmonics, Gaussian-basis 
and Lagrange-mesh methods involving different coordinate systems, are com- 
pared in studies of + bound-state energies in 3a models. Calculations are per- 
formed with different versions of the shallow Ali-Bodmer potential (with and 
| without Coulomb) and with the deep Buck-Friedrich-Wheatley potential. All 

three methods yield very accurate energies. Their advantages and drawbacks are 
evaluated. The implications of the disagreement between the obtained results 
and the experimental 12 C energies in 3a models with Coulomb interaction are 
^ ' discussed. 

> ■ 

1 Introduction 

The analysis of the structure of halo nuclei such as 6 He and n Li has renewed the 
interest for accurate calculations of three-body systems. Various dynamical approaches 
to the three-body problem are known: the hyperspherical-harmonics method Pfl |2] , the 
Faddeev coordinate-space method [Hill], the variational Gaussian-expansion method 
El HI El the Lagrange-mesh method [3 [TO] , . . . Comparing the respective merits 
of such approaches is not easy because different groups usually make different choices 
for the conditions of calculation and even for the physical constants. Moreover the 
accuracy of the presented results is not always discussed. 

Therefore, we think that it is timely to make a comparison of three of the available 
methods which are in our expertise field. Our aim is to provide results with established 
high accuracy to which other methods can be compared, under exactly the same con- 
ditions. To this end, we have selected a simple example: the 3a system with realistic 
local a-a potentials. This example contains all the difficulties of this type of problem: 
(i) Coulomb interaction between all particles, (ii) orbital-momentum dependence of the 
nuclear interaction, (iii) occurrence of forbidden states in the two-body interactions, 
(iv) full symmetrization of the three-boson system. Methods of calculation can not 
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always accurately treat all these difficulties. On physical grounds, this system is also 
interesting for several reasons. 

The three-a-boson system has a long history in nuclear physics (see Refs. fTJ ^] 
for discussions of early results). Accurate calculations could be performed with purely 
local ol-ol realistic forces of the shallow type |13j . They clearly indicated that this model 
does not provide a realistic description of the 12 C nucleus [HI HH]. Non-local two- 
body forces and/or three-body forces arising from the Pauli antisymmetrization play a 
non-negligible role. Strangely, a fair qualitative description of the bound spectrum is 
obtained when the Coulomb forces are ignored [THll5j. Therefore three-body attractive 
forces can solve the problem, at least phenomenologically ^21 0] ■ 

Deep potentials involving forbidden states are supposed to simulate the effect of 
the Pauli principle ^HlEl- Performing accurate calculations with such forces turned 
out to be much more complicated because the two-body forbidden states must be 
eliminated from the three-body calculation. Otherwise, many unphysical states show 
up [T5] . Attempts to solve that problem have led to contradictory results [THJ 1201 HE1 121] ■ 
Because of computer time limitations, whether the results obtained in Ref. [21] by one 
of us have reached convergence could not be established. It is thus timely to reexamine 
this question with several accurate methods and to establish firmly the 3a bound states 
energies. 

The 3a system has also recently gained a lot of interest as a possible example of a 
new phenomenon: the possible Bose-Einstein condensation of a bosons in the light nu- 
clei 12 C, 16 0, ... In Ref. 122], on the basis of microscopic calculations, it was proposed 
that the astrophysically significant state of the 12 C nucleus at the excitation energy 
7.65 MeV has a dominant relative s-wave structure as expected for a boson conden- 
sate. Moreover, a Generator Coordinate Method study [2H] of the 12 C and 16 nuclei 
structure predicts the existence of near-threshold na states which would be the analog 
for finite systems of an a condensation in infinite matter. However, recent variational 
calculations based on point-like bosons in 3a and 4a models with a local a-a force ob- 
tain rather small s-wave components and lead to the conclusion that an interpretation 
as a Bose-Einstein condensation is too far from reality for these nuclei [23]. In this 
physical issue, it is important to work with as exact as possible wave functions. The 
accuracy of the results must be ensured to validate any physical interpretation. 

The aim of the present paper is to compare the results and merits of three dif- 
ferent approaches of the 3a bound-state problem for two types of effective a-a in- 
teractions. The compared methods are the three-dimensional Lagrange-mesh method 
(LMM) [H1E]|, the hyperspherical-harmonics method (HHM) on a Lagrange mesh [25] . 
and the variational method on a Gaussian basis (VGM) j7[ [23- The LMM is under 
some conditions very accurate and can help calibrating the other two. The HHM is 
representative of variational techniques in hyperspherical coordinates. The VGM and 
its extension the stochastic variational method are versatile tools which can also treat 
larger numbers of particles. 

The selected effective a-a interactions are the shallow AB potential of Ali and Bod- 
mer with a strong repulsive core ^3] and the deep BFW potential of Buck, Friedrich 
and Wheatley [T7]. For elastic a-a scattering, these potentials are essentially equiv- 
alent [26J. However, not only are their form factors very different but their use for a 
three-body system lead to quite different difficulties. Since the effects of the Pauli prin- 
ciple affecting the twelve nucleons can not be treated in an exact way in a 3a model, 



one can try to simulate them with a microscopically founded potential. The deep BFW 
potential possesses unphysical bound states in the lowest s and d waves of relative mo- 
tion to simulate the Pauli-forbidden states. However these forbidden states can not be 
kept in the three-body calculation because it becomes very difficult to find the physical 
states among the many eigenvalues of the Hamiltonian ^B]- The forbidden states can 
be approximately eliminated from the solutions of the three-body Schrodinger equa- 
tion with the method of orthogonalising pseudopotentials [27j. However, this 3a model 
only takes into account two-body Pauli effects due to nucleon exchanges between two 
a particles. Using the full three-body Pauli projector is a very complicated problem 
[i?T| which will not be considered here. 

The 3a model is described in section 2, with emphasis on the potential choice. In 
section 3, the three methods of resolution of the three-body Schrodinger equation are 
summarized. In section 4, the results and the respective merits and limitations of the 
methods are compared. Some comments are made in section 5. Concluding remarks 
are presented in section 6. 



2 Model and potentials 

The Hamiltonian of three identical bosons with mass m a has the simple form 

H = "^-(V? + V 2 . + V 2 .) + V(r 23 ) + V(r 31 ) + V(r 12 ), (1) 

where 7*y = Tj — T\ is the relative coordinate between particles j and i. For deep 
two-body potentials with forbidden states in some partial waves, we use the method 
of orthogonalising pseudopotentials [27J EE] Each two-body potential V(rjj) of the 
Hamiltonian is replaced by the corresponding pseudopotential 

V(ri j ) = V(r ij )+Aj2tf j >, (2) 

f 

where the sum runs over the forbidden states. The constant A must be taken large 
enough to push the forbidden states to high energies. The projector on a two-body 
forbidden state f is given by 

ry=£l¥ , fmf)(¥'fmfl> ( 3 ) 

rrif 

where <£>f mf (r^) is the wave function of a forbidden state between a particles i and 
j. As shown below, this wave function must be accurate. Because of the symmetry 
properties of identical bosons, only even partial waves of the a-a potentials contribute. 

We use h 2 /m a = 10.4465 MeV fm 2 and e 2 = 1.44 MeV fm in all calculations. The 
Coulomb interaction is taken as 

V c {r) =4e 2 erf (f3r)/r. (4) 
The first model of the a-a interaction used in this work is the AB potential d of 

Ref. nm, 

Va.b(t) = V\ exp(-r/ir 2 ) + V 2 exp(-7? 2 r 2 ) + V c (r). (5) 



The s-wave parameters are V\ = 500 MeV, rji = 0.49 fm -2 and V2 = —130 MeV, 
i] 2 = 0.225625 fm -2 for potential d. In the simplified s-wave version (ABdo), these 
parameters are used in all partial waves and we choose (3 = 0.75 fm -1 like in the BFW 
potential below. In the /-dependent potential (ABd), the parameters are as in the s 
wave, except V x = 320 MeV for I = 2, V x = for / > 4 and (3 = y/3/(2 x 1.44) in the 
Coulomb potential for all / jT3|. The second model makes use of the BFW potential of 
Ref. UU, 

^BFw(r) = K>exp(-r /0 r 2 ) + V c {r), (6) 

with Vq = -122.6225 MeV, r/ = 0.22 fm~ 2 and (3 = 0.75 fm -1 . Both potentials describe 
fairly well the experimental phase shifts of the a-at scattering for I = 0, 2, 4 up to about 
20 MeV. The BFW potential has three nonphysical bound states forbidden by the Pauli 
principle at the energies E(0f) = -72.625691755 MeV, E(0%) = -25.618638588 MeV 
and E{2 + ) = -22.000501732 MeV. The Of state corresponds to the forbidden shell 
configuration s 8 , while the 0^ and 2 + states correspond to s 6 p 2 . 

3 Methods of resolution of the three-body Schro- 
dinger equation 

Basically, the three methods that we consider are all variational calculations of the 
three-body bound-state energies, but with different coordinate systems on one hand and 
different types of basis on the other hand. Two of the methods contain an additional 
Gauss-quadrature approximation. The three methods are based on a minimization of 
the expectation value of the Hamiltonian over linear or nonlinear variational parameters 
in a trial wave function. The comparison of their results is thus a significant test. 

The LMM |3 E5] makes use of three-dimensional Lagrange basis functions, i.e. 
infinitely differentiable functions which vanish at all points of an associated mesh, with 
the exception of one. Because of the Gauss quadrature associated with the mesh, the 
potential matrix in the Lagrange basis is diagonal and its diagonal elements are the 
potential values at mesh points. Remarkably the Lagrange-mesh method appears to 
be as accurate as the corresponding variational calculation [2H1 EE] • However, it does 
not apply to /-dependent potentials. 

The Lagrange basis functions are defined in the system of perimetric coordinates 

x = r i2 - r 32 + ri 3 , 

y = r 12 + r 32 -r 13 , (7) 
z = ~ri2 + r 3 2 + r 13 , 

where = |r» — rj\ is the distance between two particles. The basis functions read 
^ 0+ (x,y,z) =Y,C ilk Sf t (x/h)f J (y/h)f k (z/h) (8) 

ijk 

where S is the symmetrization projector and h is a scaling parameter. The one- 
dimensional Lagrange functions read 

f t (x/h) = (-l)\h Ui )^^^e-^ h (9) 



where L^{u) is the Laguerre polynomial of degree N and Ui is one of its zeros, i.e. 
Ln(ui) = 0. The interest of this approach is that the mesh equations resulting from 
the Gauss approximation are rather simple and only involve potential values at the 
(hui, huj, huk) mesh points (see Ref. P for details). 

In the following the calculations are performed with up to iV = 38 mesh points for 
each of the three dimensions. The scaling factor h is around 0.5 fm for 0^ and 1.2 fm 
for 0^" (the results are independent of the precise value of h). After symmetrization, 
the size of the largest matrix is iV(iV + 1)(N + 2)/6 = 9880 but computing times are 
short because the filling rate is about 10 %. 

In the HHM [UEIES], the variational wave function is expanded over hyperspherical 
harmonics as 

M/° + ( P ,fi 5 ) = Y,c lKi sy l L(n 5 )kp/h), (io) 

IKi 

where are hyperspherical harmonics depending on five angular variables noted 
as O5. The hyperradial wave functions depend on the rotationally and permutation- 
ally invariant hyperradius p and are expanded on regularized Lagrange-Laguerre basis 
functions 

fi{p/h) = {p/hmf^fiip/h) (11) 

where ft is given by Eq. (jHJ). Here also a Gauss approximation eliminates the need 
for calculations of potential matrix elements and significantly reduces the computing 
times (see Ref. [22] for details and tests). One of the main advantages of this method 
is its validity for both bound-state and scattering problems. 

Here we use hyperspherical harmonics up to K max = 30 and N = 30 hyperradial 
mesh points with a scaling factor h = 0.3 fm. Because of symmetrization, the basis 
size is N(K max + 2)(K mSLX + 6)/16 = 2160. 

The VGM El El EE] is a high-accuracy method for the study of the structure of 
quantum-mechanical few-body systems. Its combination with stochastic methods jHJ 
makes it applicable to few-body problems (up to six clusters at present [2H]) i n atomic, 
nuclear and quark physics. The variational wave function reads 

* 0+ (x, y) = £ C Uj S[Y l (x) ® Y l (y)) 00 x l y l exp(- ai x 2 - (3 3 y 2 ) (12) 

lij 

where x and y are the Jacobi coordinates of any of the three possible sets, and the 
on and j3j are non-linear parameters (see Ref. [21] for details and for the choice of 
non- linear parameters) . 

Matrix elements of the Hamiltonian are calculated analytically for /-independent 
potentials with simple form factors such as ABdo (see pHIE])- A numerical evaluation 
is necessary for the moments of potentials depending on I or with a complicated form 
factor. However, in most cases, one can use recurrence formulas [21] which allow to 
tabulate these moments. The high flexibility of a many-particle Gaussian basis makes 
it possible to describe several-particle configurations that are formed in the ground 
and excited states of multicluster systems [7j. For /-dependent potentials, extensions 
of the code are required with respect to Ref. [2*T] . If the potential term depends on 
a specific Jacobi coordinate x, it is necessary to perform a change of variable in the 
basis functions appearing in the matrix element and depending on another set of Jacobi 
coordinates in order to express them in the corresponding (x, y) set. After an analytic 



integration over y and over the angular components of x, the remaining integration 
over x is performed numerically. A drawback of the VGM is the non-orthogonality of 
the basis which may restrict the basis size. In the VGM calculations below, the number 
of Gaussians is 680, except otherwise indicated. 

4 Results 

The three numerical methods are used to determine energy values for the ground and 
lowest excited bound states of the 3a system with J 71 " = + . 

The numerical results obtained with the ABdo potential are presented in Table 1. 
As in Ref. let us first discuss calculations without the Coulomb term. The three 
methods yield close numbers which prove a high accuracy of all of them. The LMM is 
faster and more accurate. Its accuracy can be estimated by varying the scaling factor 
h and the number N of mesh points. ^From these tests, all displayed digits should be 
correct. The LMM accuracy is much better for the ground state than for the excited 
state. The error of the VGM is about 2 x 1(T 6 MeV. With 372 Gaussians, the results 
are -5.12205 and -1.3523 MeV while with 280 Gaussians they are -5.1215 and -1.341 
MeV. The HHM accuracy is 10~ 3 MeV for the ground state but is less good for the 
excited state. To improve this accuracy, higher values of X max should be used. 

With the Coulomb term in ABdo, only one weakly bound ground state is obtained 
(see Table 1). The convergence of all methods is slower in that case. The accuracy is 
better than 10" 8 for the LMM, 10" 5 for the VGM and 10~ 3 for the HHM, in MeV. The 
3a energy is far from the experimental 12 C energy —7.275 MeV [30J because off-shell 
effects are not well reproduced by the potential. Other Coulomb terms lead to similar 
results. With the Coulomb parameter j3 from the AB paper, one finds —0.61733886. 
For the point Coulomb potential — 4e 2 /r, the energy is —0.57238115. With the value 
h 2 /m a = 10.36675 of Ref. [21, we obtain -5.18093389 and -0.62106627 MeV without 
and with point Coulomb, respectively, in excellent agreement with the values —5.18 
and —0.62 obtained by these authors. 

Results with the full ABd potential are also presented in Table 1. The LMM does 
not apply to /-dependent potentials. We think that the accuracies of VGM and HHM 
are comparable to the ABdo case. They agree within 10 -3 MeV. The binding energy 
is increased by about 1.3 MeV (without Coulomb) or 1 MeV (with Coulomb) with 
respect to the ground-state energy of the ABd potential. In spite of this increase, the 
obtained energy remains not very realistic for 12 C. 

Energies of the 3a ground state calculated with the BFW potential including the 
Coulomb term are presented in Table 2 for the HHM and VGM. Because of the pro- 
jection on forbidden states, the LMM can not be applied. To show the convergence of 
the orthogonalising-pseudopotentials method, we present numerical results for several 
values of the projection parameter A (10 5 , 10 6 , 10 7 , 10 8 and 10 9 MeV). It was shown in 
Ref. [21 1 that lower A values can not fully eliminate the effect of the forbidden states. 
The VGM calculations are here performed with 280 Gaussians. Including more Gaus- 
sians leads to a numerical instability for large A values. However, this rather small 
Gaussian basis already yields a good estimation for the ground-state energy as shown 
by the close agreement between the two variational approaches for each A value. They 
agree within one or two percent. A fair convergence with respect to A is obtained for 



values higher than 10 7 MeV. Extrapolation with respect to K max in the HHM yields 
E = -0.30 for A = 10 7 , -0.28 for 10 8 and -0.27 for 10 9 (in MeV). The best result of 
Ref. [21 (—0.283 for A 7 = 7 and A = 10 8 ) is not far from these values. The apparent 
lack of convergence observed in that work is due to a comparison with calculations 
with an insufficient accuracy on the forbidden states. 

The 3a binding energy is 0.3 MeV smaller with the BFW potential than with the 
ABd potential and 1.25 MeV smaller than with the ABd potential. The convergence 
of our results allows us to draw a conclusion that was not accessible with the basis sizes 
in Ref. j2lj: a deep local potential accurately reproducing the a-a phase shifts does not 
provide a better description of the 3a system than shallow potentials. Independently 
of the choice of local potential, the 3a model is not valid for the shell-model-like states 
of 12 C. 

5 A few comments on the 3a model 

First we must emphasize the quality of the results obtained by Visschers and Van 
Wageningen more than 30 years ago 14]. For ABd without Coulomb, they obtain 
—6.37 MeV for the ground-state energy, only 0.05 MeV about our result. This result 
is excellent taking into account the computer limitations of that time. With Coulomb, 
their result is less good as expected for the convergence of a smaller binding energy. 

Reliable results now exist for a deep potential. They confirm that the 3a model with 
local forces is not able to reproduce the 12 C ground-state energy. The deep-potential 
energy is even above the shallow-potential one. A similar effect is observed for the 
three-nucleon system |3T] but not for the 6 He halo nucleus [23] . 

An attractive three-body force can be used to cure this problem [T""J|"""|. The effective 
three-body force employed in Ref. [?J| depends however on the hyperradius p and is 
specific for the HHM (see also Ref. j2J. This type of form factor has little physical 
meaning. 

The origin of the weakness of the 3a model should be understandable microscopi- 
cally. Microscopic 3a-cluster models "3*2*1 [3~3~j provide results much closer to experiment. 
With an effective nucleon-nucleon interaction, they are able to reproduce fairly well the 
a-a scattering and the 12 C properties, simultaneously. In particular, the 12 C ground 
state has a compact shell-model structure which can be reproduced by overlapping mi- 
croscopic a clusters. On the contrary, the 3a-boson picture with local, deep or shallow, 
a-a forces is far from reality. States with a compact shell-model structure are not sat- 
isfactorily described by the model. Moreover, the high sensitivity of the ground-state 
binding energy of the 3a system to the description of the two-body forbidden states 
[2~"""| shows the importance of a correct treatment of the Pauli principle in cluster-model 
calculations. 

The success of 12-nucleon descriptions of the 3a system has inspired a model which 
is apparently very similar to ours but which is based on a quite different philosophy. In 
this model local a-a interactions are obtained by folding an effective nucleon-nucleon 
interaction and the forbidden states are those of the microscopic a-a norm kernel. 
These forbidden states are simple oscillator states which are used in the pseudopotential 
(H3) in place of the exact bound states of the potential. This procedure leads to an 
overbinding of the 3a system |34j . in opposition with our results. 



Using oscillator forbidden states with the BFW potential would be inconsistent. 
In the simple 3a model that we have discussed, no information about the underlying 
structure of the a particle is available. The choice of the oscillator parameter in the for- 
bidden states derived from a microscopic 12-body calculation would be quite arbitrary. 
We have tested this variant by replacing in the pseudopotential the exact forbidden 
wave functions of the BFW potential by oscillator wave functions. According to the 
choice of oscillator parameter, various ground-state energies can be obtained. For each 
choice, the energy is much lower than the energy —0.26 MeV given in Table 2 since 
the bound states of the BFW potential are only partly eliminated. In principle, with 
arbitrary forbidden-state wave functions, the pseudopotential technique can provide 
any result between the energy —240.65 MeV derived without elimination of forbidden 
states and the energy —0.26 MeV corresponding to full elimination. The more-or-less 
arbitrary choice of the wave functions of the forbidden state to be eliminated becomes 
a parameter of the model with which the experimental 12 C binding energy could be 
reproduced or which can lead to overbinding as well as underbinding. 

The preceding discussion concerns local a-a potentials. The role of non-locality in 
the interaction is not yet clearly evaluated. An attempt to elucidate this problem within 
the 3a model has been performed in Ref. [33] where a non-local interaction derived from 
a resonating-group microscopic calculation fitting the a-a scattering provides a better- 
bound 12 C ground state. However, the calculation and the elimination of forbidden 
states are directly performed with the a-a T-matrix and the corresponding potential 
obtained after this elimination (which is presumably also non-local) is not available. 

6 Conclusion 

To summarize, the energies of the + ground and lowest excited states in the 3a model 
were evaluated in the framework of three variational approaches based on different 
coordinate systems. The LMM is faster and more accurate but is restricted to potentials 
which do not depend on the orbital momentum. The VGM is flexible and accurate but 
is more difficult to apply with /-dependent potentials. The non-orthogonality of the 
basis may cause problems. The slowness of convergence with respect to K makes the 
HHM slightly less accurate than the other two. However, it can easily be adapted to l- 
dependent potentials and be extended to positive energies. The elimination of forbidden 
states with pseudopotentials shows convergence with respect to the parameter A in 
contradiction with the conclusion of Ref. [21] . In order to reach convergence, the wave 
functions of the forbidden states must be very accurate. 

On the physical side, when used in a realistic way with the Coulomb interaction, 
both AB and BFW a-a potentials yield poor results for the 3a system. The 3a binding 
energies are too small compared with the 12 C experimental value. The 12 C ground-state 
structure can be described very well with microscopic 3a models [321 EB] but it can 
not be described with three-a-boson models using realistic local interactions. The role 
of nonlocality in the two-body interaction still needs to be understood. It has recently 
started to be explored in Ref. [3*3] . 
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Table 1: + energies (in MeV) of the 3a system calculated with the s-wave (ABd ) or 
/-dependent (ABd) Ali-Bodmer potentials of d type. 



state 


LMM HHM VGM 


o 2 + 


ABd without Coulomb 
-5.122093595 -5.1219 -5.1220913 
-1.3606 -1.20 -1.3566 




ABd with Coulomb {(5 = 0.75) 
-0.58427008 -0.5836 -0.584266 


o 2 + 


ABd without Coulomb 

-6.423 -6.42285 
-1.92 -1.934 


ot 


ABd with Coulomb 

-1.523 -1.523 



[21 

[22; 

[23; 
[24; 
[25; 

[26; 

[27 

[28; 
[29; 

[30; 

[31 

[32; 
[33; 
[34' 



Table 2: 3a ground-state energy calculated with the BFW potential for different values 
of the projection parameter A (in MeV). 



A 


HHM 


VGM 


10 5 


-0.638 


-0.644 


10 6 


-0.416 


-0.426 


10 7 


-0.284 


-0.288 


10 8 


-0.261 


-0.263 


10 9 


-0.259 


-0.261 



